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Recurrent microblazar activity in Cygnus X-l? 
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Abstract. Recurrent flaring events at X- and soft gamma-ray energies have been recently reported for the galactic 
black hole candidate Cygnus X-l. The observed fluxes during these transient outbursts are far higher than what is 
observed in "normal" episodes. Here we suggest that the origin of this radiation is non-thermal and produced by 
inverse Compton interactions between relativistic electrons in the jet and external photon fields, with a dominant 
contribution from the companion star held. The recurrent and relatively rapid variability could be explained by 
the precession of the jet, which results in a variable Doppler amplification. 
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1. Introduction 

Cygnus X-l is the most extensively studied black hole can- 
didate in the Galaxy. It is a very bright X-ray binary with 
a compact object of ~ 10.1 M Q and a companion 09.7 
lab star of ~ 17.8 M Q (Herrero et al. 1995), at an esti- 
mated distance of ~ 2 kpc (e.g. Gierlihski et al. 1999 and 
references therein). As in other sources of this type, the 
X-ray emission switches between soft and hard states, be- 
ing most of the time in the latter. The spectrum in both 
states can be approximately represented as the sum of a 
blackbody plus a power law with exponential cut-off (e.g. 
Poutanen et al. 1997). During the soft state the black- 
body component is dominant and the power law is steep, 
with a photon spectral index T ~ 2.8 (e.g. Frontera et 
al. 2001). During the hard state more energy is in the 
power law component, which is then harder, with pho- 
ton indices ~ 1.6 (e.g. Gierlihski et al. 1997). The black- 
body component is usually understood as emission from 
a cold, optically thick accretion disk, whereas the power 
law component is thought to be originated in an optically 
thin hot corona by thermal Comptonization of disk pho- 
tons (Poutanen et al. 1997, Dove et al. 1997). The hot 
corona would fill the inner few tens of gravitational radii 
around the black hole. The accretion disk penetrates only 
marginally in the coronal region. In the hard state the 
thermal X-ray emission is dominated by the corona, with 
typical luminosities of a few times 10 37 erg s _1 . In the soft 
states, the disk approaches to the black hole and then most 
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of the energy dissipation occurs through it (Poutanen et 
al. 1997; see also Poutanen & Coppi 1998). Cygnus X-l has 
a persistent, mildly variable, compact continuum coun- 
terpart of flat spectrum (e.g. Pooley et al. 1999). During 
many years, evidence for non-thermal radio jets in Cygnus 
X-l was lacking, despite the efforts of the observers (e.g. 
Marti et al. 1996). Finally, the jet was detected by Stirling 
et al. (2001) at milliarcsecond resolution using VLBA ob- 
servations. The jet-like feature extends up to ~ 15 mas 
with an opening angle of less than 2 degrees. The spec- 
turn seems to be flat, and no counterjet is observed. The 
total radio emission at 8.4 GHz is ~ 11 mJy, with varia- 
tions of ~ 2 mJy over timescales of 2 days (Stirling et al. 
2001). The average angle with the line of sight, if the jet 
is perpendicular to the disk, seems to be ~ 30° (Fender 
2001). Very recently, the Interplanetary Network detected 
a transient soft-gamma ray event from the general direc- 
tion of Cygnus X-l (Golenetskii et al. 2002). Analysis of 
previous data indicates that at least other two events were 
observed during 1995. These latter events were also de- 
tected by BATSE instrument on the Compton Gamma- 
Ray Observatory, suggesting that they were originated in 
Cygnus X-l (Schmidt 2002). The luminosities above 15 
keV of the outbursts were in the range 1 — 2 x 10 38 erg 
s _1 , much higher than the typical thermal luminosity in 
the hard state. In this letter we suggest that these flaring 
events can be interpreted in terms of non-thermal microb- 
lazar activity (Mirabel & Rodriguez 1999). We study the 
effects of the interaction of the relativistic jet with the 
ambient photon fields from the accretion disk, the corona, 
and the companion star, and we calculate the expected 
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Fig. 1. Sketch of the general situation discussed in the pa- 
per. A relativistic jet is injected close to the black hole in 
Cygnus X-l. This jet must traverse photon fields created 
by the cold accretion disk, the hot corona, and the stel- 
lar companion. Inverse Compton up-scattering of some of 
these photons is unavoidable. 
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Fig. 2. External photon fields to which the jet is exposed. 



non-thermal contribution to the keV-MeV spectrum. The 
recurrent character of the events can be explained through 
variable Doppler boosting originated in the precession of 
the jet (Kaufman-Bernado et al. 2002). In the next section 
we present the model, and then we discuss the implica- 
tions. 



2. Non-thermal high-energy emission for Cygnus 
X-l 

We shall consider the effects of the injection of a rela- 
tivistic leptonic jet at a few Schwarzschild radii from the 
central black hole and its subsequent propagation through 
the ambient photon fields (see a sketch of the situation 
in Figure 1). The individual electrons have Lorentz fac- 
tors 7 in the lab frame and the flow is assumed to have 
a bulk Lorentz factor T. In accordance to the disk/jet 
symbiosis model (e.g. Falcke & Biermann 1999, Markoff 
et al. 2001) and the estimated accretion rate of Cygnus 
X-l (~ 10~ 8 M Q yr" 1 , e.g. Poutanen et al. 1997), we 
adopt a mean particle density n ~ 10 14 cm~ 3 for the 
jet at 10r s from the black hole. The electron energy dis- 
tribution will be assumed to be a power law given by 
(in the lab frame, Georganopoulos et al. 2001): n(j) = 
■^D 2+p ^ p P(j 1 D 1 72-D, 7), where k is a constant and 
D = [r (1 — /3cos0)] _1 is the usual Doppler factor: (3 and 
<f> are the bulk velocity in units of c and the viewing angle, 
respectively. The function P is 1 for 71 _D < 7 < 72 1? and 
otherwise. 

The jet will traverse first the coronal region. The parame- 
ters that characterize this region and its photon field in the 
hard state change with time. We shall assume here the typ- 
ical values given by Poutanen et al. (1997) in their Table 



1. All photon fields that interact with the jet are shown 
in Figure 2. The photon field of the corona is assumed to 
be isotropic. External inverse Compton interactions with 
these photons can be treated using the head-on approxi- 
mation in the Klein-Nishina regime since 760 ^ 1, with 
eo the average energy of the coronal photons in units of 
the electron rest energy (eom e c 2 ~ 100 keV). The reader 
is referred to Georganopoulos et al. (2001) for details of 
calculation. In the case of the photons from the accretion 
disk and the companion star we can work in the Thomson 
regime, taken into account that the photons from the disc 
come from behind the jet, introducing an additional fac- 
tor (1 — cos 0)( p+1 )/ 2 that reduces the effects of beaming 
(Dermer et al. 1992). The stellar photons, instead, can be 
treated as an isotropic field. The emerging spectrum of the 
specific luminosity can be approximated by a power law of 
index (p — l)/2 in the Thomson regime (Georganopoulos 
et al. 2001, 2002; Kaufman-Bernado et al. 2002). When 
72 3> 71 this can be written as: 

-(p-l)/2 



L f = 



dL 
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kVa T c U2P- 



(1) 



dedVL 7reo(l +p)(3 + p) \e 

where or is the Thomson cross section, V is the inter- 
action volume, and U is the average energy density of 
the photon field. In the Klein-Nishina regime, where nu- 
merical integrations are necessary, the results significantly 
depart from a power law, resulting in a softer spectrum. 
The Compton losses in the different regions will modify 
the injected electron spectrum, introducing a break in 
the power law at the energy at which the cooling time 
equals the escape time. This will occur at (e.g. Longair 
1997): 7b = (3TO c c 2 )/(4er T r 2 i7£ esc ). Here i osc is average 
time spent by the particles in the field region (typically 
^csc ~ l/c, with I the linear size), and the rest of the sym- 
bols have their usual meanings. The spectrum will steepen 
from an index p to p+ 1 for energies higher than 7b. After 
the interaction of the jet with the coronal region, the mod- 
ified spectrum will be injected in the stellar photon field 
region, suffering further losses and modifications. In our 
calculations we have adopted two different values for the 
original electron energy index: a hard index p — 1.5 and a 
steeper index p = 2.3. The jet was assumed forming an an- 
gle of 30 degrees with the line of sight, and a bulk Lorentz 
factor r = 5 was adopted. The initial part of the jet was 
modeled as a cylindrical structure; the coronal region, fol- 
lowing Poutanen et al. (1997), was considered as a spher- 
ical region of ~ 500 km in radius. Gamma-rays produced 
within this region will be mostly absorbed in the local 
field through pair creation. Using the simple formulae by 
Herterich (1974) along with the adopted parameters for 
the coronal region, we estimate an optical depth r <~ 1 for 
photons of 5 MeV. The probability for a 10 MeV photon to 
escape from the corona is only 0.1. The results obtained 
for the different models are shown in Figures 3 and 4. 
There, we show the different components due to external 
Compton scattering of the different photon fields in the 
lab frame, for the different injected electron distributions 
we have assumed. This non-thermal contribution should 
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Fig. 3. Results of the model for an injection electron spec- 
trum with index p = 1.5 in a cylindrical jet forming a 
viewing angle of 30 degrees. The bulk Lorentz factor is 
r = 5 and the electron power law extends from 71 = 2 to 
72 = 10 3 . Three different components are shown, resulting 
from the up-scattering of star (A), disk (B), and corona 
(C) photons. Radiation absorbed in the local photon fields 
is shown in dashed lines. Notice that the contribution from 
the coronal photons is not a power law because of the 
Klein-Nishina effects. 
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Fig. 5. Variation of the amplification factor for the ex- 
ternal inverse Compton emission from the jet (electron 
power-law indices of p = 1.5 -dots- and p = 2.3 -solid 
line-) as a function of time in the precessing microblazar 
model for Cygnus X-l (precessing opening angle of 16.5 
degrees). Time units are normalized to the precessing pe- 
riod. 



tected by the Interplanetary Network. The transit through 
the peak of flux magnification can be very fast in the ob- 
server's frame, leading to quick and transient states when 
the total flux is dominated by the non-thermal contribu- 
tion. The angular velocity of the tidally induced precession 
can be approximated by uj p « — 3Gmcos^/4r 3 1 cjd (Katz 
1980, Romero et al. 2000), where 9 is the half-opening an- 
gle of the precession cone, m is the mass of star, r m is the 
binary separation, and Wd is the Keplerian velocity of the 
accretion disc, which can be written in terms of the mass 
M of the black hole as w d = (GAf/r 3 ) 1 / 2 . For 9 = 16.5 
deg as we have assumed, the observed orbital period of 5.6 
days, and a precessing period of ~ 140 days (Brocksopp et 
al. 1999, see below), we get a disk size of ~ 3.9 x 10 11 cm, 
quite reasonable for a wind-accreting system like Cygnus 
X-l. 



be added to the thermal components shown in Figure 2 in 
order to recover the total emission. We can see that the 
non-thermal emission is dominated by the up-scattering of 
the stellar photons. Above 5 MeV almost all energy flux 
comes from this component. For a viewing angle of 30 de- 
grees, we still have the thermal emission dominating by a 
factor of 5. But if we introduce the gravitational effects 
of the companion star, then precession of the disk should 
occur (Larwood 1998, Kaufman-Bernado et al. 2002). In 
Figure 5 we show the modification of the beaming ampli- 
fication factor of the external inverse Compton emission 
for a precessing angle of 16.5 degrees. The time axis is 
normalized in units of the period T. We see that there is 
a variation of about 1 order of magnitude in the emission 
measured in the observer's frame because of the preces- 
sion. This means that when the jet is closer to the line 
of sight, the non-thermal luminosity can reach values of 
~ 10 38 erg s^ 1 , as observed in the recurrent outbursts de- 



3. Discussion 

The above outlined model incorporates the different 
known components of Cygnus X-l, accretion disk, corona, 
stellar companion, and relativistic jet, in an integrated pic- 
ture where transient non-thermal outbursts are a natural 
and expected result. The amplitude of these outbursts can 
be similar to what has been recently observed in some in- 
triguing flaring episodes in this source. Our model is differ- 
ent from the model recently proposed by Georganopoulos 
et al. (2002) not only because we incorporate the effects 
of precession, but also because we do not attempt to ex- 
plain the bulk of X-ray emission as non-thermal all the 
time. This emission is normally dominated by thermal 
Comptonization in the hot corona around the black hole, 
except during the microblazar phase, and in this case we 
incorporate the effects of the interaction of the jet with 
the corona in our calculations. We emphasize that, as it 
is shown in Figures 3 and 4, during the transient microb- 
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lazar phase the X- and soft 7-ray spectrum will be softer 
than in the normal hard state, when the coronal emis- 
sion dominates. This is an unavoidable consequence of the 
steepening produced by Compton losses in the injected 
electron spectrum and can be used to test our proposal, 
not only through new observations of Cygnus X-l, but 
also of other potential microblazars as LS5039 (Paredes 
et al. 2000). Recently, Brocksopp et al. (1999) have found 
multiwavelength evidence for the presence of a 142.0 ±7.1 
days period in Cygnus X-l. The optical and X-ray period 
seem to be originated in the precession of the accretion 
disc (Brocksopp et al. 1999), whereas the modulation at 
radio wavelengths is probably produced by the associated 
precession of the jet (see Pooley et al. 1999) . The morphol- 
ogy of the extended radio jet, with a clear bend, is also 
consistent with a precession of the inner beam (Stirling et 
al. 2001). The periodic signal in the radio lightcurve, how- 
ever, is not expected to be as strong as at high energies 
since the magnification factor for the synchrotron emis- 
sion goes as £>( 3 +p)/ 2 , whereas for the external Compton 
component it goes as D 2+p (Georganopoulos et al. 2002). 
The time lag between the two high-energy flares observed 
in 1995 is ~ 75 days, about a half of the value reported by 
Brocksopp et al. (1999), but since Cygnus X-l is a wind- 
accreting system variations in the period along a span of 
several years are possible. Certainly, more observations on 
longer time spans are necessary to constrain the dynam- 
ical models. In the model presented here the duty cycle 
of the blazar phase is rather small, ~ 10%. Future X-ray 
observations of non-thermal flares can be used for a better 
determination of the geometric parameters. An interesting 
feature of our model is that most of the gamma-rays pro- 
duced within the coronal region will be absorbed by pair 
production. Sooner or later these pairs will annihilate pro- 
ducing a broad, blucshiftcd feature in the MeV spectrum. 
Details of calculation of the emerging spectrum are be- 
yond the scope of this Letter (see Abraham et al. 2001), 
but it is clear that the forthcoming INTEGRAL satellite 
will be able to probe Cygnus X-l spectrum and its tem- 
poral evolution at this energy range, helping to test and 
constrain the model here proposed. Hopefully, very soon 
we will be able to clarify the role played by non-thermal 
processes in this fascinating object. 
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